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Electrochemical  impedance  analysis  is  performed  to  predict  the  rate  capability  of  two  commercial 
activated  carbon  electrodes  (RP20  and  MSP20)  for  electric  double-layer  capacitor.  To  this  end,  ac 
impedance  data  are  fitted  with  an  equivalent  circuit  that  comprises  ohmic  resistance  and  impedance  of 
intra-particle  pores.  To  characterize  the  latter,  ionic  accessibility  into  intra-particle  pores  is  profiled  by 
using  the  fitted  impedance  parameters,  and  the  profiles  are  transformed  into  utilizable  capacitance  plots 
as  a  function  of  charge-discharge  rate.  The  rate  capability  that  is  predicted  from  the  impedance  analysis 
is  well-matched  with  that  observed  from  a  charge-discharge  rate  test.  It  is  found  that  rate  capability  is 
determined  by  ionic  accessibility  as  well  as  ohmic  voltage  drop.  A  lower  value  in  ionic  accessibility  for 
MSP20  is  attributed  to  smaller  pore  diameter,  longer  length,  and  higher  degree  of  complexity  in  pore 
structure. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Electric  double-layer  capacitors  (EDLCs)  have  been  used  as  a 
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cycle  life  (>100,000  cycles).  At  present,  activated  carbons  are  the 
most  popular  EDLC  electrodes  due  to  their  high  electrical  conduc¬ 
tivity,  large  surface  area  and  wide  pore  size  distribution  [1  . 

In  general,  the  rate  capability  of  activated  carbon  electrodes  is 
governed  by  two  factors:  1)  ohmic  resistance  that  determines  the 
ohmic  voltage  drop  and  thus  working  voltage,  and  2)  ionic  acces¬ 
sibility  into  pores,  by  which  utilizable  capacitance  is  determined. 
Commonly,  dc  methods  are  employed  to  assess  the  rate  capability; 
for  instance,  galvanostatic  charge-discharge  cycling  with  varied 
current  density.  Here,  the  ohmic  resistance  (R0 hm)  can  be  estimated 
from  the  voltage  drop  at  the  beginning  of  current  reversal,  whereas 
the  ion  accessibility  can  be  estimated  from  the  delivered  capaci¬ 
tance  values.  Even  if  the  dc  methods  can  give  general  information 
on  EDLC  parameters,  ac  methods  such  as  electrochemical  imped¬ 
ance  spectroscopy  (EIS)  can  provide  more  detailed  information;  for 
instance,  pore  structure  of  electrode  materials  and  kinetics  in 
double-layer  charging/discharging  processes  [2  .  In  the  previous 
works,  ac  impedance  analysis  based  on  transmission-line  model 
with  pore  size  distribution  (TLM-PSD)  was  successfully  utilized  to 
assess  the  pore  structure  of  carbon  materials  [3-5  .  Later  on,  the 
concept  of  TLM-PSD  has  been  incorporated  into  complex  capaci¬ 
tance  analysis,  which  allows  a  graphical  analysis  on  capacitance, 
rate  capability,  and  leakage  current  [6-15  .  Especially,  it  was 
demonstrated  that  the  effect  of  pore  structure,  electrode  potential, 
and  electrode  thickness  on  rate  capability  can  be  easily  evaluated 
by  comparing  the  peak  frequency  on  imaginary  capacitance  plots. 

The  primary  objective  of  this  work  is  to  establish  an  analytical 
method  to  correlate  the  ac  method  (impedance  analysis)  with  dc 
method  (galvanostatic  charge-discharge  cycling)  for  the  prediction 
of  rate  capability  of  EDLC  systems.  To  this  end,  experimental  ac 
impedance  data  are  analysed  to  obtain  ohmic  resistance  and  ionic 
accessibility  for  activated  carbon  electrodes.  Then,  the  utilizable 
capacitance  is  calculated  as  a  function  of  operating  time  (top)  or 
current  density  ( i )  by  using  correlation  functions  between  ac  and  dc 
variables.  Finally,  the  transform  technique  is  validated  by 
comparing  the  rate  capability  predicted  from  two  methods. 

2.  Experimental  details 

2.1.  Material  characterizations 

The  samples  were  characterized  by  using  a  field-emission 
scanning  electron  microscope  (FE-SEM,  JEOL  JSM-6700F),  nitro¬ 
gen  adsorption  (Micromeritics,  ASAP  2010),  and  small  angle  X-ray 
scattering  (SAXS,  Bruker  GADDS,  CuKa,  A  =  0.154056  nm).  The  ni¬ 
trogen  adsorption  data  was  analysed  by  Barret-Joyner-Halenda 
(BJH)  and  modified  micropore  (MP)  16,17  methods  to  characterize 
mesopores  and  micropores,  respectively. 

2.2.  Electrode  preparation 

For  the  electrochemical  tests,  the  composite  electrodes  were 
prepared  by  coating  the  slurry  of  active  material  (RP20  or  MSP20), 
polytetrafluoroethylene  and  carboxyl  methyl  cellulose 
(PTFE  +  CMC,  6:4  in  mass  ratio)  binder,  and  conductive  carbon 
(Super-P)  (8:1 :1  in  mass  ratio  in  deionized  water)  on  a  piece  of  Al  foil 
(thickness  =  21  pm).  The  electrode  plates  were  dried  in  vacuum 
oven  at  120  °C  for  12  h  without  pressing  process.  The  resultant 
thickness  of  the  coated  film  was  approximately  40  ±  2  pm  with  the 
active  mass  loading  of  1.5  ±  0.1  mg  on  the  electrode  area  of  0.95  cm2. 

2.3.  Electrochemical  measurements  and  impedance  fitting 

For  a  3-electrode  cell  test,  a  home-made  cell  comprising  poly¬ 
ether-ether-ketone  (PEEK)  body  and  stainless  steel  (316L)  current 


collectors  was  used.  A  spring  was  attached  to  a  current  collector  to 
keep  the  positive,  negative,  and  reference  electrodes  be  tightly 
contacted  to  the  current  collectors.  O-rings  (Viton®)  were  used  to 
ensure  the  sealing  at  the  joints  of  current  collectors  and  cell  body. 
Positive  and  negative  electrodes  were  positioned  to  exactly 
confront  each  other,  and  the  reference  electrode  was  just  beside 
them.  Two  pieces  of  separator  was  put  between  positive,  reference, 
and  negative  electrodes  for  all  the  electrodes  to  be  insulated  one 
another.  Activated  carbon  (RP20)  were  utilized  to  prepare  a  refer¬ 
ence  electrode  18]  (coated  film  on  Al  foil,  1  mg  on  0.6  cm2,  50  pm 
thickness)  and  a  counter  electrode  (sheet-type,  30  mg  on  1.8  cm2, 
400  pm  thickness).  A  porous  glass  fibre  was  used  as  the  separator. 
The  electrolytes  were  1  M  tetraethylammonium  tetrafluoroborate 
(TEABF4)  in  acetonitrile  (AN)  or  propylene  carbonate  (PC).  Imped¬ 
ance  was  measured  at  0.0  V  vs.  carbon  over  the  frequency  range  of 
2  mHz-100  kHz  (Zahner,  Im6e)  with  a  root  mean  square  (rms) 
amplitude  of  5  mV.  For  the  complex  nonlinear  least  squares  (CNLS) 
fitting  of  the  impedance  data,  the  TLM-PSD  was  coded  in  FORTRAN 
to  run  in  LEVM  8.09  software  [19  .  Modulus  weighting  (with 
respect  to  calculated  values)  was  adopted  for  the  CNLS  fitting. 

For  rate  experiment,  symmetric  2-electrode  (CR2032  coin  type) 
cells  were  assembled  by  sandwiching  two  identical  electrodes  and 
operated  in  0-3.5  V.  Experimental  conditions  are  same  with  the  3- 
electrode  experiments,  except  for  the  electrode  area  (2.27  cm2). 
Galvanostatic  charge-discharge  cycling  was  made  with  a  WBCS- 
3000  battery  cycler  (Wonatech  Co.).  Both  of  charging  and  dis¬ 
charging  current  densities  were  varied  from  0.5  to  40  mA  cm-2. 

3.  Results  and  discussion 

3.2.  Characterization  of  activated  carbons 

As  the  EDLC  electrodes,  two  different  activated  carbons  were 
used;  RP20  (Kuraray  Chemical  Co.)  and  MSP20  (Kansai  Coke  and 
Chemicals  Co.).  The  former  is  produced  by  physical  steam  activa¬ 
tion,  whereas  the  latter  by  chemical  activation  using  potassium 
hydroxide  (KOH).  As  shown  in  the  FE-SEM  images  (inset  of  Fig.  la), 
RP20  has  a  rougher  morphology  as  compared  with  that  for  MSP20. 
It  is  known  that  steam  activation  (RP20)  leads  to  more  severe 
morphological  change  as  compared  with  KOH  activation  (MSP20) 
[20].  The  particle  diameter  is  similar  for  two  carbons  (~5  pm). 

From  the  nitrogen  adsorption  isotherms,  total  pore  volumes 
(VWai)  are  calculated  to  be  0.77  cm3  g-1  (RP20)  and  0.98  cm3  g-1 
(MSP20)  ( Fable  1).  The  adsorption-desorption  isotherms  (Fig.  la) 
also  show  that  RP20  has  a  larger  portion  of  mesopores  and  macro¬ 
pores,  while  the  micropores  are  dominant  in  MSP20.  The  BET  sur¬ 
face  area  (SBet)  is  larger  for  MSP20  due  to  the  higher  population  of 
micropores.  In  the  BJH  pore  size  distribution  (Fig.  lb),  RP20  shows  a 
larger  portion  of  meso-  or  macropore  volume  between  the  di¬ 
ameters  of  4-300  nm.  The  meso-  and  macro-pore  volume  at 
D  >2  nm  are  0.17  cm3  g-1  for  RP20  and  0.15  cm3  g-1  for  MSP20, 
which  is  22%  and  15%  of  Vtotai,  respectively.  Accordingly,  the  average 
mesopore  diameter  (DmeSo)  is  calculated  to  be  larger  for  RP20 
(4.9  nm)  as  compared  with  that  for  MSP20  (2.8  nm).  When  the 
micropore  region  (D  <  2  nm)  is  analysed  by  the  modified  MP  method 
(Fig.  lc),  larger  portion  of  micropore  volume  is  calculated  for  MSP20 
(0.83  cm3  g-1,  85%)  compared  to  RP20  (0.60  cm3  g_1,  78%). 

Pore  complexity  is  compared  for  two  activated  carbons  by  using 
SAXS  technique  (Fig.  Id).  In  the  Porod  region,  the  slope  in  loga¬ 
rithmic  scale  (Porod  slope,  z)  is  related  with  the  fractal  dimension 
(d)  of  pores  by:  Intensity  (q)  ~  q-(6-d)  =  q~z.  The  theoretical  Porod 
slope  is  4  for  an  ideal  flat  surface,  but  the  slope  gradually  decreases 
to  approach  1  with  an  increase  in  pore  complexity  [21].  From  the 
Porod  slopes  in  this  work  (ca.  2.2  for  RP20  and  1.4  for  MSP20),  it  can 
be  assumed  that  MSP20  has  a  more  complex  pore  structure. 
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Fig.l.  (a)  N2  adsorption  isotherms  and  pore  size  distribution  by  (b)  BJH  and  (c)  modified  MP  methods  for  RP20  and  MSP20.  (d)  SAXS  patterns  of  RP20  and  MSP20.  FE-SEM  images  of 
RP20  and  MSP20  in  the  inset. 


Based  on  the  above  results,  it  is  expected  that  RP20  exhibits 
higher  rate  capability  since  it  has  larger  pore  size  and  lower  pore 
complexity  22  ,  but  lower  specific  capacitance  as  its  surface  area  is 
smaller  than  that  for  MSP20.  Meanwhile,  if  the  specific  capacitance 
is  calculated  with  an  assumption  that  the  whole  surface  is  utilized 
for  charge  storage  with  identical  capacitance  of  8  pF  cm-2  [23,24],  it 
is  estimated  to  be  120  F  g-1  (RP20)  and  160  F  g_1  (MSP20)  from  the 
BET  surface  areas.  However,  the  calculated  value  is  only  a  rough 
estimation  since  the  pore  utilization  and  unit  capacitance  are 
strongly  dependent  on  the  used  electrolyte  and  charge/discharge 
rate  in  practical  EDLCs. 

32.  Graphical  analysis  of  EIS  data:  Nyquist  plot  and  complex 
capacitance  analysis 

Fig.  2  displays  the  Nyquist  plots  that  are  obtained  from  RP20  and 
MSP20  in  1.0  M  TEABF4/AN  and  TEABF4/PC  electrolyte,  in  which  a 
semicircle  and  a  spike  appear  at  the  high  and  low  frequency  region, 
respectively.  Total  impedance  can  be  divided  into  three  compo¬ 
nents  in  the  Nyquist  plots;  bulk  electrolyte  resistance  (x-intercept 
at  the  highest  frequency  region),  interfacial  impedance  between 
electrode  and  bulk  solution  (semicircle  at  the  middle  frequency 
region),  and  the  impedance  that  is  associated  with  intra-particle 
pores  (spike  at  the  low  frequency  region)  (Fig.  3).  The  first  two 
terms  are  mainly  dependent  on  the  electrolyte  solution,  while  the 
last  one  (low  frequency  tails)  is  controlled  by  both  electrode  ma¬ 
terials  and  electrolytes. 


Table  1 

Pore  properties  derived  from  nitrogen  adsorption  isotherms. 


SBET/m2g  la 

VWai/cm3  g  1 

DaVg/nm 

f^meso/nrn 

RP20 

1520  ±  30 

0.77 

2.0 

4.9 

MSP20 

2060  ±  40 

0.98 

1.9 

2.8 

a  Measured  by  Brunaur— Emmett— Teller  (BET)  method. 
b  Measured  by  Barret— Joyner— Halenda  (BJH)  method. 


Previously,  the  semicircle  at  the  middle  frequency  region  has 
been  assigned  to  the  ion  transport  processes  25-27]  or  the  contact 
impedance  between  electrode  and  current  collector  28,29].  As  the 
resistances  are  larger  in  the  PC-based  electrolyte  for  both  elec¬ 
trodes  in  this  work,  the  ion  transport  processes  seem  to  be  domi¬ 
nant  in  our  experiment,  even  if  the  contribution  from  the  electronic 
contact  resistance  cannot  be  neglected.  The  semicircle  has  also 
been  assigned  to  the  inter-particle  ionic  impedance  [25  ,  but  our 
simulation  shows  that  the  inter-particle  ionic  resistances  are 
merged  in  the  sloped  line  at  low  frequencies  as  additional  pore 
resistances  [9  .  Practically,  the  process  for  the  semi-circle  can  be 
regarded  as  a  simple  resistance  at  lower  frequency  than  10  Hz, 
while  the  typical  EDLC  charge/discharge  condition  corresponds  to 
the  frequency  of  0.1  Hz. 


Fig.  2.  Nyquist  plots  of  (a)  RP20  and  (b)  MSP20  electrodes  in  AN  or  PC  electrolytes. 
Lines:  CNLS  fitted  results. 
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Fig.  3.  Equivalent  circuit  model  for  the  porous  carbon  electrodes. 


As  a  simple  approach,  capacitive  electrochemical  systems  can  be 
represented  by  a  serial  connection  of  equivalent  series  resistance 
(Re sr)  and  utilizable  capacitance  (Cutii),  which  approximates  the 
capacitor  behaviour  at  sufficiently  low  frequencies.  In  this 
approximation,  Res r  comprises  1)  bulk  electrolyte  resistance,  2) 
interfacial  resistance,  and  3)  apparent  resistance  of  intra-particle 
pores  (Fig.  3).  The  capacitive  charge  storage  is  assumed  to  occur 
on  electrochemically  equivalent  sites  with  identical  serial  resis¬ 
tance.  As  the  impedance  of  such  circuit  is  Res r  -  j/(wCutu),  where  j  is 
imaginary  unit  and  w  is  angular  frequency  (=2tu/),  the  Res r  and  Cutii 
can  be  determined  from  the  real  and  imaginary  part  of  experi¬ 
mental  impedance  data  from  the  perpendicular  region  in  the 
Nyquist  plots  (Fig.  2). 

In  the  AN-based  electrolyte  (ionic  conductivity  =  56  mS  cm-1  at 
25  °C),  from  the  impedance  data  at  10  mHz,  the  Res r  and  Cutii  are 
calculated  to  be  21  Q  cm2  and  80  F  g-1  for  RP20,  and  17  Q  cm2  and 
140  F  g-1  for  MSP20.  In  the  PC-based  electrolyte  (ionic 
conductivity  =  13  mS  cm'1  at  25  °C),  the  Resr  and  Cutji  are  calculated 


to  be  42  Q  cm2  and  75  F  g'1  for  RP20,  and  47  Q  cm2  and  130  F  g-1  for 
MSP20.  Here,  mass  (g)  refers  to  mass  of  carbon  only  ( i.e .  without  the 
mass  of  aluminium  current  collector).  As  higher  rate  capability  can 
be  assumed  by  smaller  time  constant  (r)  that  is  the  product  of  Re sr 
and  Cutii,  it  is  predicted  that  the  rate  capability  of  RP20  (r:  2.6  s  in 
AN  and  4.7  s  in  PC)  is  higher  as  compared  with  that  for  MSP20  (r: 
3.5  s  in  AN  and  9.1  s  in  PC).  Rate  capability  can  be  predicted  from 
this  simple  graphical  analysis  by  using  Resr  and  Cutii  values,  but 
detailed  analysis,  such  as  the  separate  discussion  on  the  effect  of 
ion  size  and  conductivity,  is  not  available. 

More  detailed  analysis  on  the  capacitive  parameters  and  rate 
capability  can  be  performed  by  using  the  complex  capacitance 
analysis  [7].  For  this,  the  measured  impedance  (Z (J)  =  Z'(f)  +jZ"(f )) 
is  transformed  into  the  complex  capacitance  (C (J)  =  C'(f)  +  jC"(f))  by 
the  relation  of  C(f)  =  1  /[/wZ(f)].  The  real  Part  °f  complex  capacitance 
(C'C/))  represents  the  apparent  capacitance  value  as  a  function  of 
frequency,  while  the  imaginary  part  (C'(/))  is  correlated  with  C'(f)  by 
Kronig-Kramers  (K-K)  relations.  When  C"(f)  is  plotted  as  a  func¬ 
tion  of  frequency  in  semi-log  scale,  peak-shaped  curves  are 
appeared  in  the  imaginary  capacitance  plot  (C"(f)  vs.  log/,  Fig.  4a 
and  b).  Then,  as  suggested  by  the  K-I<  relations,  total  capacitance 
(Ctot)  can  be  calculated  from  the  peak  area  (Ap)  as  Ctot  =  1.466AP 
[13].  In  addition,  the  peak  frequency  (/p)  at  the  maximum  C"(f) 
corresponds  to  the  characteristic  frequency,  which  is  inversely 
proportional  to  the  time  constant  (r)  of  capacitive  systems. 
Therefore,  capacitance  and  rate  capability  can  be  easily  estimated 
from  this  graphical  analysis. 

Fig.  4a  and  b  presents  the  imaginary  capacitance  plots  obtained 
from  two  activated  carbon  electrodes.  The  Ctot  and  /p  values  are 
determined  to  be  83.2  F  g'1/70-7  mHz  (RP20)  and  140.4  F  g-1/ 
50.4  mHz  (MSP20)  in  the  AN-based  electrolyte;  and  79.2  F  g-1/ 
34.9  mHz  (RP20),  and  143.0  F  g_1/16.6  mHz  (MSP20)  in  the  PC- 
based  electrolyte  ( Table  2).  The  Ctot  values  are  comparable  in  two 
electrolytes,  but  they  were  smaller  than  the  estimations  from  the 
BET  surface  area  (i.e.  120  and  160  F  g  1  for  RP20  and  MSP20, 
respectively;  assuming  8  pF  cm-2  for  carbon  surface)  for  both  RP20 
(-69%)  and  MSP20  (-89%)  because  very  small  pores  that  are  not 
utilizable  for  ion  adsorption  are  also  detected  by  N2  molecules 
(BET).  This  finding  supports  that  EDLC  can  be  more  properly  eval¬ 
uated  by  electrochemical  method  that  utilizes  actual  ions  as  the 
probe  molecule,  compared  to  N2  adsorption  method  that  uses  N2  as 
the  probe  molecule  [5].  Meanwhile,  the  peak  frequency  values 
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Fig.  4.  Imaginary  capacitance  plots  of  (a)  RP20  and  (b)  MSP20  electrodes  in  AN  or  PC  electrolytes  with  CNLS  fitted  results.  Ionic  accessibility  profiles  for  (c)  RP20  and  (d)  MSP20 
electrodes  with  an  assumption  of  log-normal  distribution  (Eq.  (6)). 
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Table  2 

Capacitive  parameters  derived  from  imaginary  capacitance  plots. 


/p/mHz 

Ctot/F  g-1 

FWHM 

RP20 

1  M  TEABF4/AN 

71 

83 

1.33 

1  M  TEABF4/PC 

35 

79 

1.36 

MSP20 

1  M  TEABF4/AN 

50 

140 

1.33 

1  M  TEABF4/PC 

17 

143 

1.35 

indicate  that  rate  capability  of  RP20  is  higher  than  that  of  MSP20  in 
both  electrolytes  by  a  factor  of  1.4  (AN)  and  2.1  (PC).  This  solvent- 
dependent  rate  capability  for  two  carbon  electrodes  can  be 
ascribed  to  the  difference  in  the  solvated  ionic  diameter:  ca.  1.2  nm 
for  AN  and  1.4  nm  for  PC  [30  . 


CtlM— PSd(/)  =  Qot  J  C°(/,  a0)p(a0)d\n  a0  (5) 

—  oo 

Based  on  Eq.  (5),  the  ionic  accessibility  profile,  p(a0),  can  be 
obtained  by  de-convoluting  the  impedance  data  with  discrete 
Fourier  transform  [7,33].  However,  it  has  been  reported  that  the 
PSD  of  porous  materials  can  be  assumed  to  have  log-normal  func¬ 
tion  (Eq.  (6))  to  reasonably  fit  the  experimental  EIS  data,  where  a0* 
and  a  represent  the  characteristic  ionic  accessibility  and  the  degree 
of  distribution  for  the  electrode,  respectively  [3-5  . 


P(a  0) 


3.3.  Ionic  accessibility  profile  by  CNLS  fitting 

A  schematic  Nyquist  plot  is  presented  in  Fig.  3a.  At  the  high 
frequency  limit,  a  resistive  term  appears  due  to  ion  transport  in 
bulk  solution,  which  can  be  represented  by  a  simple  resistance  of 
ftbuik  (Fig-  3b).  The  semicircle  in  the  middle  frequency  region  is 
assigned  to  the  ion  migration  across  the  bulk-electrolyte/electrode 
interface,  while  the  sloping  spike  at  the  low  frequency  region  is  for 
the  combined  effect  of  ion  transport  and  double-layer  formation 
inside  pores.  As  these  processes  are  serially  connected,  total 
impedance  (Z .(f))  is  the  sum  of  bulk  resistance  (Rbuik),  impedance 
for  the  semicircle  (Zinterface)  and  sloping  spike  (Zintra-Pore): 

Z(/)  =  ^bulk  +  ^interface  (f)  "P  ^intra-pore(f)  0  ) 

To  simulate  the  interfacial  impedance  (Zinterface),  (AiQiX^O) 
circuit  (Fig.  3c)  is  used,  where  Qis  the  constant  phase  element  (CPE, 
Zcpe  =  l/[F(/w)p])  [31  .  The  R2C2  is  added  because  the  depressed 
semicircles  are  not  successfully  fitted  by  the  R1Q1  alone.  In  the 
fitting,  p  is  ca.  0.9  for  all  the  fitted  results  and  T  is  a  variable.  For  the 
Zintra-pore,  the  cylindrical  pore  model  is  utilized,  in  which  the 
equivalent  circuit  follows  the  TLM  with  segmental  ionic  resistances 
and  surface  capacitances  in  intra-particle  pores  [32  .  For  a  single 
pore  or  uniform  multiple  pores,  the  electrochemical  characteristics 
with  ac  signals  (CTlm(/,^o))  can  be  expressed  as  the  product  of  total 
capacitance  (Ctot)  and  the  characteristic  function  (C °(f,a0))  as: 

CTLM(/,«o)  =  CtotxC0(f,«o)  (2) 

C°(f,«o)=-^=tanh(^)  (3) 

vW  V  “°  / 

a°~2  VQP  (4) 

Here,  the  ionic  accessibility  (a0)  indicates  the  feasibility  for  ac 
signal  to  penetrate  into  pores  (radius:  r,  and  length:  l)  when  the 
ionic  conductivity  in  pores  is  k  and  the  capacitance  per  unit  area  on 
pore  surface  is  Cd  (Eq.  (4)).  For  a  pore  with  an  a0,  the  frequency 
dependency  of  the  capacitance  utilization  can  be  represented  by 
Eq.  (3). 

The  uniform  pore  model  fails  to  analyse  practical  porous  car¬ 
bon  electrodes  since  pore  structures  are  significantly  non- 

uniform,  except  for  some  highly  ordered  materials  [3  .  There¬ 

fore,  it  is  necessary  to  consider  the  non-uniformity  of  pores  in 
modelling  the  electrochemical  characteristics  of  porous  carbon 
electrodes.  When  p(a0)  is  the  ionic  accessibility  profile  as  a  func¬ 
tion  of  ionic  accessibility,  the  complex  capacitance  of  non-uniform 
multiple  pores  (CTLM-psdC0),  with  total  capacitance  of  Ctot,  is  given 
as  follows  [7]: 


The  experimental  impedance  data  of  RP20  and  MSP20  elec¬ 
trodes  are  CNLS  fitted  by  Eq.  (1),  assuming  non-uniform  multiple 
pores  with  log-normal  distribution  (Eq.  (6)).  The  Zintra-pore(fl  is 
represented  as  1/L/wCtlm-psdC/)]  (Eq.  (5)).  The  fitted  curves  are  in  a 
good  agreement  with  the  experimental  data  as  shown  in  Nyquist 
plots  (Fig.  2)  and  imaginary  capacitance  plots  (Fig.  4a  and  b).  From 
the  optimized  parameters  ( fable  3),  the  electrochemical  charac¬ 
teristics  for  intra-particle  pores  can  be  separately  analysed  from 
that  of  the  bulk  electrolyte  and  interface  (Table  4). 

When  the  fitted  parameters  are  compared  with  those  from  the 
complex  capacitance  analysis,  total  capacitance  values  are  similar 
within  1%  of  difference.  From  the  larger  a0*  values,  higher  rate 
capability  is  expected  for  RP20  in  both  electrolytes,  which  is  in 
accordance  with  the  expectation  from  the/p  values  in  the  graphical 
analysis.  Therefore,  it  can  be  concluded  that  the  capacitance  and 
rate  capability,  which  are  the  key  parameters  to  characterize  elec¬ 
trochemical  characteristics  of  porous  carbon  electrodes,  can  be 
successfully  analysed  by  graphical  analysis  with  C"(f)  vs.  log  /plot, 
as  confirmed  by  fitted  parameters. 

The  ionic  accessibility  profile,  which  enables  the  quantification 
of  rate  capability,  is  derived  by  the  unique  interpretation  of  the 
impedance  data.  The  ionic  accessibility  profile  is  defined  as  the 
distribution  density  of  capacitance  existence  with  respect  to  the 
natural  log  of  ionic  accessibility  (p(a0)  vs.  In  a0),  which  is  obtained 
from  Eq.  (6)  with  the  fitted  parameters  a0*  and  a  (Fig.  4c  and  d).  The 
distribution  function,  p(a0)  =  d[C/Ctot]/d[ln  a0],  can  be  interpreted 
as  the  distribution  densities  of  surface  existence  (d[S/Stot]/d[ln  a0]) 
if  constant  Cd  is  assumed  over  the  entire  surface.  As  the  ionic 
accessibility  profile  of  RP20  is  shifted  to  the  higher  a0  direction,  it 
can  be  noticed  that  RP20  has  more  accessible  pores  by  ac  signals 
compared  with  MSP20.  In  other  words,  MSP20  has  smaller  a0*, 
which  is  in  accordance  with  the  nitrogen  adsorption  (larger  portion 
of  micropores)  and  SAXS  analysis  (higher  complexity  of  porous 
structure)  results.  The  smaller  pore  diameter,  more  complex  pore 
structure,  and  possibly  longer  average  pore  length  of  MSP20  result 
in  the  smaller  a0*  value.  Of  two  electrolytes,  the  AN-based  one 
provides  larger  a0*  due  to  higher  ionic  conductivity  as  compared 
with  the  PC-based  one.  All  the  properties  are  colligated  to  a0*,  the 
characteristic  ionic  accessibility  into  pores,  according  to  Eq.  (4).  The 
larger  values  of  a  for  RP20  reflect  that  the  pore  structure  is  more 
largely  distributed  compared  with  MSP20. 


Table  3 

TLM-PSD  parameters  obtained  by  fitting  the  impedance  data  to  Eq.  1 


Ctot/Fg-1 

<c/s-03 

a 

x2 

RP20 

ANa 

82.1  ±  0.2 

1.21  ±  0.09 

1.22  ±  0.07 

7.2  x  10~3 

PC 

79.0  ±  0.1 

0.68  ±  0.01 

1.15  ±  0.03 

9.1  x  10~3 

MSP20 

AN 

139.0  ±  0.2 

0.56  ±  0.02 

0.73  ±  0.03 

7.9  x  10~3 

PC 

141.5  ±  0.2 

0.273  ±  0.001 

0.622  ±  0.009 

4.2  x  10  3 

a  All  the  electrolytes  contain  1  M  TEABF4. 
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Table  4 

Parameters  corresponding  to  the  impedance  of  bulk  electrolyte  and  interfaces  obtained  by  fitting  the  impedance  data  to  Eq.  1 


Rbuik/Q  cm2 

Ri/Q  cm2b 

R2/Q  cm2 

Td  10~6 

Pi 

C2/(iF  cm  2 

RP20 

ANa 

0.37  ±  0.01 

9.7  ±  0.4 

4.7  ±  0.3 

68  ±  6 

0.89  ±  0.01 

420  ±  60 

PC 

1.04  ±  0.03 

20.6  ±  0.6 

6.9  ±  0.6 

37  ±  2 

0.89  ±  0.01 

260  ±  30 

MSP20 

AN 

0.43  ±  0.02 

4.4  ±  0.5 

5.1  ±  0.5 

70  ±  10 

0.88  ±  0.02 

80  ±20 

PC 

1.02  ±  0.02 

19.6  ±  0.6 

9.5  ±  0.6 

46  ±  2 

0.88  ±  0.01 

140  ±  20 

a  All  the  electrolytes  contain  1  M  TEABF4. 

b  Rn  stands  for  interfacial  resistance  of  nth  semicircle,  7i  and  pi  for  CPE  parameters  of  1st  semicircle,  and  C2  for  interfacial  capacitance  of  2nd  semicircle. 


3.4.  Differential  and  apparent  capacitances  calculated  from  ionic 
accessibility  profiles 

From  the  ionic  accessibility  profiles  (Fig.  4c  and  d),  the  rate 
capability  of  two  activated  carbon  electrodes  can  be  predicted.  To 
this  end,  firstly,  the  utilization  of  cylindrical  pores  with  an  ionic 
accessibility  of  a0,  at  an  operating  ac  frequency  of/op,  is  assumed  to 
be  Re[C°(/0p^o)]  based  on  the  TLM  (Eq.  (3)).  Then,  with  a  decrease 
in  /0p,  the  utilization  of  a  single  pore  will  gradually  increase  from 
0  (high  frequency  limit)  to  1.0  (low  frequency  limit).  Theoretically, 
the  pore  utilization  will  be  10%,  50%,  and  90%  at/op/a2  of  15.9, 0.991, 
and  0.294,  respectively  (Fig.  SI). 

As  the  porous  carbon  electrodes  have  non-uniform  pores,  the 
differential  specific  capacitance  of  an  electrode,  as  a  function  of/op, 
can  be  calculated  from  the  EIS  data  (Cdiff.EisC/op))  according  to  the 
following  integral  equation  (Eq.  (7)),  where  the  p(a0)  represents  the 
ionic  accessibility  profile  obtained  from  EIS  analysis,  and  Ctot  the 
average  total  capacitance  of  an  electrode  in  the  operating  voltage 
range. 


Re 


C0  (/op>  OLq 


p(a0) din  a0 


If  the  capacitance  is  constant  over  the  potential  range,  the  Ct0 1 
can  be  determined  from  the  Ctot  value  measured  by  single  EIS 
analysis  at  any  potential.  The  most  direct  method  for  studying  the 
potential-dependent  quantities  (e.g.  capacitance)  is  multiple 
impedance  measurement  at  every  potential,  which  is  not  practical 
because  each  impedance  measurement  takes  at  least  2  h  for  a 
particular  potential  point.  Alternatively,  either  cyclic  voltammetry 
(CV)  at  a  slow  scan  rate  or  charge-discharge  cycling  at  a  slow  rate 
can  be  performed  to  measure  Ct0 t.  In  this  study,  using  slowly 
scanned  CV,  Ctot  values  are  determined  to  be  120  F  g  1  for  RP20  (in 
AN  and  PC)  and  160  F  g_1  and  150  F  g_1  for  MSP20/AN  and  MSP20/ 
PC,  respectively  (Fig.  5).  As  the  capacitance  value  is  potential- 
dependent,  the  measured  Ctot  Are  1.1 -1.5  times  larger  than  the 
Qot  value  obtained  from  EIS  analysis  at  0  V  vs.  carbon.  Usually  the 
capacitance  of  activated  carbon  electrodes  is  minimum  at  the  po¬ 
tential  of  zero  charge  (pzc,  about  0  V  vs.  carbon)  [34].  The  lower 
specific  capacitance  of  MSP20  in  PC  compared  to  AN  is  due  to  the 
shrinkage  of  the  voltammogram  at  about  <-1.5  V  vs.  carbon,  which 
has  been  explained  by  either  pore-narrowing  by  ion-induced  strain 
or  saturation  of  the  pores  by  ions  before  the  electrode  potential 
reaches  the  lower  limit  [35,36  . 

Fig.  6a  and  b  shows  the  Cdiff,Eis(top)  curves  as  a  function  of 
operating  time  (top)  for  RP20  and  MSP20  electrodes  in  AN-  and  PC- 
based  electrolytes.  The  operating  frequency  (/op)  is  transformed  into 
the  operating  time  (top)  in  dc  techniques  by  top  =  0.5 //op,  considering 
the  period  of  ac  signal  as  the  time  span  for  a  charge  and  discharge 
cycle.  The  Cdiff,Eis(£op)  curve  predicts  the  differential  capacitance  as  a 
function  of  top.  For  example,  the  differential  capacitance  of  an  EDLC 
cell  with  RP20/AN  is  expected  to  be  120  F  g-1  at  slow  charge/ 
discharge  (top:  1000  s)  and  decreased  to  109  F  g-1  at  higher  rate 
(10  s).  In  the  case  of  MSP20  with  PC  electrolytes,  the  decrease  in 


differential  capacitance  can  be  predicted  to  be  more  severe: 
150  F  g-1  (top:  1000  s)  to  93  F  g_1  (10  s).  Accordingly,  the  “rate 
capability”  can  be  quantitatively  determined  for  various  combina¬ 
tions  of  electrode  materials  and  electrolytes.  When  the  operation 
time  of  the  EDLC  cells  are  decreased  from  1000  s  to  10  s,  the  retained 
capacitances  will  be:  RP20/AN  (90.8%)  >  MSP20/AN  (85.8%)  >  RP20/ 
PC  (82.7%)  >  MSP20/PC  (62.0%).  In  general,  the  higher  the  a0*  value, 
the  more  feasible  is  the  operation  at  shorter  top.  In  addition  to  the 
effect  of  a0*,  the  Cdiff,Eis(top)  curves  change  less  steeper  for  RP20 
with  respect  to  top,  due  to  the  larger  a  values  (1.2)  than  those  of 
MSP20  (0.6-0.7).  Also,  the  smaller  value  of  a  for  MSP20  leads  to  the 
higher  susceptibility  to  the  ionic  conductivity  for  the  rate  capability, 
which  is  evidenced  by  the  larger  difference  of  Cdiff,Eis(top)  curves  by 
electrolytes  for  MSP20  (0.64  decade)  compared  to  that  of  RP20  (0.50 
decade)  (Fig.  6a  and  b). 

To  validate  the  developed  EIS  analysis  method  in  a  voltage  range 
as  wide  as  possible,  symmetrical  EDLC  cells  were  constructed  and 
galvanostatic  charge-discharge  measurements  were  carried 
out  between  0  and  3.5  V  at  various  current  densities 
(i  =  0.5-40  mA  cm-2)  (Fig.  7).  The  EDLC  cells  were  evaluated  at  an 
operating  voltage  of  3.5  V,  which  was  used  by  some  groups  for 
higher  energy  density  [37,38  .  Note  that  the  3.5  V  operation  was 
used  for  proof  of  concept,  and  the  practical  operating  voltage  of 
EDLC  is  between  2.5  and  2.8  V,  which  can  provide  better  durability. 

The  differential  specific  capacitances  for  an  electrode,  Cdiff,di- 
sch(0>  were  determined  from  the  discharge  capacity  of  a  symmetric 
full  cell  (Qce n,disch(0)  by  using  following  equation.  Here,  the  Cdiff,di- 
sch(0  is  inversely  proportional  to  the  slope  in  the  linear  region  that 
appears  after  the  IR- drop.  In  Eq.  (8),  Qceii.dischO)  was  multiplied  by  4 


E  /  V  vs  carbon 


E  /  V  vs  carbon 


Fig.  5.  Cyclic  voltammograms  of  (a)  RP20  and  (b)  MSP20  electrodes  in  AN  and  PC 
electrolytes  (v  =  1  mV  s”1). 
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Fig.  6.  (a,b)  Cdiff, eis(  top)  profiles  (lines)  with  respect  to  the  operational  time  (top)  and  Cdiff.dischttop)  values  from  the  full-cell  rate  test  (points).  (c,d)  Qeii.EisO)  profiles  (dashed  lines, 
0  =  1;  solid  lines,  0  =  0.492)  with  respect  to  the  current  density  (i)  and  the  full-cell  rate  test  data  (Cce n.dischO),  points). 


considering  that  the  symmetric  EDLC  cell  contains  active  mass  of 
both  electrodes  (i.e.  total  mass  becomes  2  times  of  an  electrode's 
mass)  while  the  cell's  capacitance  becomes  half  of  an  electrode's 
capacitance  due  to  the  serial  connection. 

Qliff,disch(0  —  4  x  Qce ll,disch (0 / (^op  —  ^ohm)  (8) 

As  shown  in  Fig.  6a  and  b,  the  Cdiff, disch(  fop)  values,  which  were 
measured  by  the  galvanostatic  charge-discharge  of  EDLC  cells, 
were  well  matched  with  the  Cdiff, eis(  fop)  curves  predicted  by  the  EIS 
data  in  half-cell  tests.  Here,  the  operational  time  is  equal  to  the 
experimental  discharging  time.  For  example,  the  Cdiff  of  an  EDLC 
electrode  with  RP20/AN  was  experimentally  measured  to  be 
119  F  g-1  (0.5  mA  cm-2),  110  (10  mA  cm-2),  and  100  (40  mA  cm-2), 
which  well  agrees  with  the  expected  values  by  EIS  analysis.  For  four 
kinds  of  EDLCs,  the  deviation  between  experimental  and  predicted 
values  was  within  10%.  Therefore,  it  can  be  concluded  that  the  rate 
capability  of  differential  capacitance  of  EDLC  cells  can  be  quanti¬ 
tatively  predicted  by  the  developed  analysis  technique  of  the 
impedance  data  for  porous  electrodes  in  half  cell  test. 


Fig.  7.  Charge-discharge  curves  at  various  current  density  (0.5-40  mA  cm  2)  for  (a) 
RP20  and  (b)  MSP20  electrodes  in  AN  or  PC  electrolytes. 


In  addition  to  the  differential  specific  capacitance  of  an  elec¬ 
trode  (Cdiff),  the  apparent  specific  capacitance  of  EDLC  cells  (Cceii) 
and  its  dependency  on  the  current  density  are  another  important 
factor  in  evaluating  practical  EDLC  cells.  The  Cceii  represents  the 
practical  utilizable  capacitance  of  EDLC  cells  that  includes  the  effect 
of  ohmic  losses,  whereas  the  Cdiff  indicates  the  characteristics  of 
electrode  materials.  When  charge/discharged  at  various  currents 
(i),  the  apparent  specific  capacitance  of  an  EDLC  cell  (Cceii,disch)  is 
calculated  by  dividing  the  discharging  specific  capacity  of  sym¬ 
metric  full-cell,  Qceii,disch(i)»  by  the  operating  voltage  (Vop,  3.5  V  in 
this  work)  as: 

Q:ell, disch  (0  =  Qcell,disch(0/^op  (9) 

At  slow  current  density  (5  mA  cm-2),  measured  Cceii,disch(0 
values  for  EDLC  cells  (points  in  Fig.  6c  and  d)  are  close  to  one 
fourth  of  the  Cdiff, disch  W  values  for  single  electrodes,  as  expected 
with  negligible  i'R-drop.  With  current  increase,  the  Cceii, disch  values 
are  decreased  by  increased  ohmic  drop  (Vohm)  as  well  as  by  the 
decrease  in  Cdiff,disch-  As  a  result,  the  decrease  with  current  is  ex¬ 
pected  to  be  larger  for  Cceii,disch  than  Cdiff, disch-  For  example,  when 
the  discharge  current  was  increased  from  0.5  mA  cm-2  to 
40  mA  cm-2,  the  Cceii, disch  was  decreased  by  74%,  while  the 
decrease  in  Cdiff, disch  was  42%.  As  the  iR- drop  plays  an  important 
role,  the  effect  of  electrolytes  on  rate  capability  was  much  larger 
for  the  apparent  capacitances  ( Cceii,disch( 0 )-  When  the  electrolyte 
was  changed  from  AN  to  PC,  the  Cceii, disch  of  RP20  was  decreased 
by  19%  (40  mA  cm-2),  while  the  corresponding  Cdiff, disch  change 
was  5%. 

The  apparent  specific  capacitance  of  an  EDLC  cell,  Cceii,Eis,  also 
can  be  predicted  quantitatively  from  the  Cdiff, eis( fop)  of  an  electrode 
determined  through  the  EIS  analysis  (Eq.  (10)).  In  Eq.  (10),  Cdiff, ei- 
s(fop)  was  divided  by  4  considering  that  the  symmetric  EDLC  cell 
contains  2  times  of  active  mass  (i.e.  mass  of  activated  carbon)  while 
the  capacitance  itself  becomes  half  due  to  the  serial  connection.  For 
this,  the  ohmic  resistance  of  EDLCs  and  operating  time  should  be 
determined  at  various  current  densities. 


r  (\\  _ rt  (\\  l\r  _ Oiiff,EIS (^op)  vop 

Cceii, EIS  w  —  C'cell,EIs(I)/V/op  — - ~r - X 


Vr 


^ohm 


V( 


(10) 


op 
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The  ohmic  resistances  of  EDLC  cells  are  mainly  composed  of  bulk 
resistance  in  organic  electrolytes  and  resistive  terms  from  two 
electrodes.  In  the  EIS  analysis  with  a  three  electrode  cell,  the 
interfacial  impedance  (Zinterface)  of  a  porous  carbon  electrode  was 
represented  by  overlapped  semicircles.  However,  as  active  fre¬ 
quency  of  the  semicircle  (100  kHz-10  Hz)  is  sufficiently  larger  than 
the  operational  condition  (i.e.  current  density)  of  EDLCs 
(/op  <  0.5  Hz),  the  interfacial  impedance  (Zinterface)  can  be  regarded 
as  a  simple  resistance,  interface-  Therefore,  the  ohmic  resistance  of 
an  EDLC  cell  (R0hm)  can  be  calculated  from  the  EIS  data  as  the  sum 
of  ftbuik  and  interfacial  resistances  of  both  electrodes 
(^interface  =  2(Ri  +  R2)),  with  an  empirical  correction  factor  (0),  as: 


^ohm  —  [^bulk  +  2(i?i  +  ^2)]  x  0  (H) 

If  the  fabrication  condition  of  EDLC  cells  is  identical  to  that  of 
half-cell  test,  0  will  become  unity.  However,  the  effective  bulk  re¬ 
sistances  and  interfacial  resistances  can  be  influenced  by  the 
electrode  fabrication  methods  and  cell  geometry  [39  .  In  such  cases, 
the  correction  factor  0  is  expected  to  be  determined  empirically 
and  can  be  utilized  in  the  prediction  of  the  EDLC  performance  from 
EIS  data  according  to  the  analysis  technique  developed  in  this 
study. 

For  each  combination  of  porous  carbons  and  electrolytes,  the 
Rohm  value  was  determined  by  plotting  the  ohmic  voltage  drops  in 
the  charge-discharge  data  as  a  function  of  current  density  (Fig.  S2). 
Then,  a  linear  correlation  between  R0 hm  values  of  symmetrical 
EDLCs  and  EIS  analysis  results  (Rbuik  +  2(Ri  +  R2))  were  found 
among  MSP20/AN,  RP20/AN,  and  RP20/PC,  and  the  correction  fac¬ 
tor  0  for  the  used  EDLC  cells  was  determined  to  be  0.492  (Fig.  S3). 
This  result  implies  that  the  applied  pressure  was  probably  higher 
for  the  coin-type  cell  (EDLCs  for  charge/discharge)  compared  to  the 
home-made  test  cell  (half-cells  for  EIS  analysis).  Even  though  the 
applied  pressure  was  not  numerically  determined,  it  can  be 
concluded  that  the  effect  of  cell  fabrication  method  on  the  ohmic 
resistance  values  was  maintained  similarly  in  this  study  and 
therefore  the  prediction  of  R0hm  from  impedance  data  is  possible 
with  a  predetermined  correlation  factor. 

In  the  case  of  MSP20/PC,  the  voltage  drop  at  the  charge/ 
discharge  reversal  and  resultantly  calculated  R0hm  was  much  larger 
than  that  expected  from  the  EIS  result  with  0  =  0.492.  This  seems  to 
be  originated  from  the  additional  voltage  drops  at  the  negative 
electrode,  as  experimentally  confirmed  only  for  the  charge/ 
discharge  reversal  of  MSP20/PC  (Fig.  8).  Previously,  similar  phe¬ 
nomenon  was  reported  for  nanoporous  activated  carbon  electrode 
in  TEABF4/PC  electrolyte,  which  was  explained  by  strong  adhesion 
of  cations  on  nanopores  [40].  Even  though  the  effect  of  abnormal 
adhesion  cannot  be  fully  predicted  by  EIS  analysis,  such  combina¬ 
tions  of  porous  carbon  and  electrolytes  with  abnormal  ohmic  drop 
will  be  not  suitable  for  practical  applications  as  EDLC  devices  with 
high  efficiency. 

The  operating  time  (top)  that  corresponds  to  the  charge/ 
discharge  at  various  current  densities  (i)  could  be  determined  by 
utilizing  Eq.  (12).  When  top  is  large  (slow  charge/discharge  with  low 
i),  it  will  be  inversely  proportional  to  the  current  density  (i),  where 
the  voltage  drop  by  ohmic  drop  is  negligible  and  Cdiff, eis(  fop)  is 
constantly  maintained.  In  contrast,  with  fast  charge/discharge  (high 
i),  top  decline  with  higher  i  became  more  rapid  due  to  both  the 
significant  ohmic  drop  and  decrease  in  Cdiff, eis(  fop)-  The  relationship 
between  top  and  i  are  numerically  calculated  from  Eq.  (12b)  and 
presented  in  Fig.  S4,  where  Uop  is  3.5  V  and  R0hm  is  determined  by 
Eq.  (11 )  with  </>  =  0.492.  It  is  noted  that  at  i  <  10  mA  cm”2,  the  top  vs.  i 
profiles  of  identical  activated  carbon  is  very  similar,  regardless  of 
the  used  electrolytes. 
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Fig.  8.  Potential  profiles  of  both  electrodes  during  the  galvanostatic  charge-discharge 
of  full  cell  ( i  =  0.5  mA  cm  2).  Grey:  cell  voltage.  Red  and  blue:  potential  of  (+)  and  (-) 
electrode,  respectively.  (For  interpretation  of  the  references  to  colour  in  this  figure 
legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


n  rs\  _Cdiff,EIs(M  ( 

(2cell,EIs(l)  - - 7 - x  (/op 


A 

i^ohm)  —  i  x  fop  x  —  (12a) 


^op 

^ohm  3"  4topm  xAj Cdiff  ,EIS  (fop) 


(12b) 


Here,  A  is  the  apparent  electrode  area  (e.g.  2.27  cm2  in  this  rate 
experiment)  and  m  is  the  total  active  mass  in  a  symmetric  EDLC. 

Fig.  6c  and  d  compared  the  apparent  specific  capacitances 
calculated  from  EIS  data  in  comparison  to  the  experimental  values. 
Solid  lines  in  Fig.  6c  and  d  represent  thus  calculated  Cceii,Eis(i) 
profiles  after  calibrating  the  R0hm  values  with  the  coefficient  of 
0  =  0.492.  The  Cceii,Eis(0  and  Cceii,  disch(0  are  matched  within  6% 
deviation  for  RP20/AN,  RP20/PC,  and  MSP20/AN.  When  the  effect  of 
cell  fabrication  and  geometry  was  not  considered  (dashed  lines, 
0  =  1)  the  calculated  values  became  more  deviated  from  the 
experimental  values,  which  confirms  that  the  appropriate  correc¬ 
tion  factor  is  required  to  utilize  EIS  data  for  EDLC  cells.  The  MSP20/ 
PC  cell  shows  larger  deviation  (up  to  33%)  even  if  the  R0hm  value  was 
calibrated  with  0  =  0.492,  which  can  be  accounted  for  by  limitation 
in  estimating  the  voltage  drop,  as  described  above. 

Additionally,  charge-discharge  curves  can  be  directly  simulated 
from  the  Cdiff, eis(Q  profiles  and  calibrated  R0hm  values  at  a  current 
density  (Fig.  9),  which  are  in  agreement  with  experimentally  ob¬ 
tained  curves  except  for  MSP20/PC.  Note  that  the  simulated  curves 
were  solely  obtained  by  EIS  parameters  and  Ctot  parameter  from 
cyclic  voltammetry.  The  simulated  and  experimental  voltage  pro¬ 
files  are  matched  within  3%  deviation  for  RP20  cells,  and  4%  devi¬ 
ation  for  MSP20/AN  cell.  Due  to  the  abnormal  adhesion  of  PC- 
solvated  TEA+  ion  in  MSP20,  much  larger  deviation  (41%)  was 
observed  for  the  MSP20/PC  cell. 


3.5.  Applications  of  EIS  analysis  of  EDLCs  for  the  rate  capability 
prediction  and  design  of  capacitor  cells 

The  rate  capability  from  Cdiff, eis(  fop)  is  expected  to  be  the 
intrinsic  properties  of  activated  carbon  particles  as  the  extrinsic 
properties  are  separated  by  fitting  analysis.  In  other  words,  this 
information  is  independent  of  the  extrinsic  parameters  of  elec¬ 
trodes  (i.e.  due  to  the  electrode  preparation  process),  as  they  are 
derived  solely  from  the  intrinsic  parameters  of  intra-particle  pores; 
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Fig.  9.  Comparison  of  experimental  and  simulated  (f  =  0.492)  galvanostatic  voltage 
profiles  at  i  =  20  mA  cm  2. 

the  precedent  impedances  that  are  subject  to  the  experimental 
conditions  are  factored  out  by  fitting  analysis.  And  it  will  be 
possible  to  distinguish  intrinsic  and  extrinsic  effects  on  the  net  rate 
capabilities  of  cells  that  are  fabricated  by  different  methods. 

The  CdifffEis(top)  or  Cceii,Eis(top)  plots  can  be  applied  to  estimate 
the  rate  capability  of  a  specific  cell  or  to  design  cells  according  to 
the  desired  specification  ( i.e .  operational  rate  and  capacitance).  If  an 
EDLC  cell  is  made  of  MSP20  (10  mg)|l  M  TEABF4  in  AN  |  RP20 
(15  mg),  the  capacitance  of  positive  (C+(top))  and  negative  (C_(top)) 
electrodes  can  be  calculated  by  multiplying  active  mass  to  the 
Qiff.Eis(fop)  profiles  in  Fig.  6a  and  b.  And  the  capacitance  of  the  full¬ 
cell  (Cfc(toP))  can  be  calculated  according  to  the  rate  conditions,  by 
using  l/Cfc(top)  =  1  /C_|_( top )  +  1/C  -(fop)- 

On  the  other  hand,  the  Cdiff.Eis(fop)  curve  can  be  utilized  to  design 
a  cell  with  desired  rate  capability.  When  1  Ah  of  capacity  is  required 
at  40  mA  cm-2  for  a  3  V  operating  symmetric  full-cell  with 
R0hm  =  15  Q  cm2,  total  active  mass  of  porous  carbons  can  be  esti¬ 
mated  to  be  55  g  (RP20)  or  41  g  (MSP20)  to  fabricate  the  cell  with 
AN  electrolyte  from  Eq.  (8).  If  R0hm  is  reduced  to  2  Q  cm2,  the 
required  mass  becomes  45  g  (RP20)  or  33  g  (MSP20).  This  calcu¬ 
lation  is  also  possible  for  the  PC-based  electrolyte  at  various 
discharge  rates  utilizing  the  Cdiff.Eis(top)  plots.  It  is  noteworthy  that, 
while  we  try  to  balance  the  active  mass,  the  rate  capability  can  be 
also  lowered  with  thickness  increase  by  several  tens  of  pm  9], 
where  the  calculated  active  mass  by  the  impedance  analysis  should 
be  regarded  as  the  minimum  mass  of  active  materials  under  no 
thickness  effect.  Related  future  work  is  to  develop  more  sophisti¬ 
cated  model  that  considers  thickness  effect. 

As  this  paper  aims  to  represent  fundamental  ideas  for  estima¬ 
tion  of  rate  capability,  the  experiments  are  simplified.  Even  if  the 
principles  will  be  the  same  as  described  in  this  paper,  there  are 
some  specific  details  to  consider  before  practical  application.  (1) 
The  rate  capability  depends  on  electrode  thickness  as  well  as 
electrode  material  [9,41-44  ,  but  the  former  effect  is  not  consid¬ 
ered  in  this  study.  As  the  effect  of  electrode  thickness  can  be  ana¬ 
lysed  by  similar  impedance  analysis  [9],  more  precise  cell  design 
will  be  possible  through  this  approach. 

(2)  Also  it  is  noteworthy  that  the  ionic  accessibility  profiles  are 
obtained  from  the  impedance  measurements  at  pzc,  and  the 
possible  differences  between  the  positive  and  negative  electrodes 
are  not  represented.  If  positive  and  negative  electrodes  are  to  be 
studied  separately,  impedance  should  be  measured  at  the  middle  of 
the  operational  potential  range  for  positive  and  negative  electrodes, 
respectively.  This  approach  can  be  a  future  extension  of  this  work. 


4.  Conclusions 

Through  the  impedance  analysis,  rate  capabilities  of  porous 
carbon  electrodes  are  quantitatively  profiled  as  ionic  accessibility 
profile,  which  is  described  as  the  distribution  of  capacitance  with 
respect  to  ionic  accessibility  (a0).  The  ionic  accessibility  profiles  are 
transformed  into  calculated  specific  capacitance  of  an  electrode 
(Qiff,Eis(fop))  with  respect  to  operational  time  (top).  More  practically, 
calculated  apparent  specific  capacitance  of  a  cell  (Cceii,Eis(0)  was 
derived  from  the  calculated  specific  capacitance  of  an  electrode 
(Oiiff,Eis(top))  and  calibrated  R0hm-  The  calculated  specific  capaci¬ 
tance  of  an  electrode  (Cdiff,Eis(£op))  and  of  a  cell  (Cceii,Eis(0)  are  in  a 
good  agreement  with  the  galvanostatic  full-cell  rate  test  results. 
This  confirms  that  the  rate  capabilities  of  porous  carbons  can  be 
quantitatively  predicted  by  EIS  analysis,  considering  both  intrinsic 
(e.g.  Cdiff,Eis(top))  and  extrinsic  (e.g.  R0hm)  factors.  The  calculated 
specific  capacitance  (Cdiff,Eis(£op))  profile  can  be  used  for  assessing 
and  designing  EDLC  cells. 

The  intrinsic  and  extrinsic  factors  that  control  the  EDLC  char¬ 
acteristics  are  separated  by  their  characteristic  frequencies  or 
phase-shifts,  and  obtained  quantitatively  by  fitting  the  impedance 
data  to  appropriate  equivalent  circuits.  This  separation-ability  of 
the  impedance  analysis  may  lead  to  the  standardized  evaluation  of 
porous  carbons  for  EDLC  electrodes. 

The  ionic  accessibilities  of  porous  carbons  are  represented  by 
characteristic  ionic  accessibility  (a0*)  and  the  degree  of  distribution 
(er).  The  pore  diameter,  existence  of  subnano-pores  or  large  meso- 
pores,  fractal  dimension,  and  average  pore  length  are  the  factors 
that  control  the  ct0*.  The  smaller  the  a  value,  the  higher  is  the 
susceptibility  to  the  electrolyte  properties  (e.g.  ionic  conductivity) 
for  the  rate  capability.  For  an  electrode  with  small  a  (e.g  MSP20 
electrode),  proper  choice  of  electrolyte  is  indispensable  for  the  rate 
capability. 
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Glossary 

A:  apparent  electrode  area  (cm2) 

AN:  acetonitrile 
Ap:  peak  area 

a0:  penetrability  coefficient  (s-0  5) 

a0*:  characteristic  penetrability  coefficient  (s-0  5) 

BFJ:  tetrafluoroborate  anion 
C:  complex  capacitance  (=C  +  jC") 

C':  real  part  of  complex  capacitance  (Re[C]) 

C":  imaginary  part  of  complex  capacitance  (Im[C]) 

Cdiff,Eis-  differential  specific  capacitance  of  an  electrode  by  EIS  analysis  (F  g-1) 

Q eii.Eis-  apparent  specific  capacitance  of  a  symmetric  full  cell  by  EIS  analysis  (F  g-1) 
CNLS:  complex  nonlinear  least  squares 
CPE:  constant  phase  element 

CjLMdoio):  complex  capacitance  of  a  single  TLM  element 

Ctlm-psd ■'  complex  capacitance  of  non-uniform  multiple  pores  described  by  TLM-PSD 
C°(fa0):  characteristic  function  of  a  single  TLM  element 
Ctot:  total  capacitance  of  an  electrode  (F  g_1) 

Cm  •  the  average  total  capacitance  of  an  electrode  in  the  operating  voltage  range  ( F  g“ 1 ) 
Cutn:  utilizable  capacitance  (F  g-1) 

D;  pore  diameter  (nm) 

EIS:  electrochemical  impedance  spectroscopy 
fop:  operational  frequency  (Hz) 
fp:  peak  frequency  (Hz) 

<fi:  correction  factor  of  ohmic  resistance 

i:  current  density  (mA  cm-2) 

j:  imaginary  unit  (=>/^T) 

m:  total  active  mass  in  a  symmetric  EDLC 

p(a0):  ionic  accessibility  profile 

PC:  propylene  carbonate 

pzc:  potential  of  zero  charge 

Rn:  interfacial  resistance  of  nth  semicircle 

Cn:  interfacial  capacitance  of  nth  semicircle 

Qn.-  CPE  of  nth  semicircle 

Rbuik •'  bulk  resistance  (Q  cm2) 

Resr:  equivalent  series  resistance  (Q  cm2) 

Rinterface:  interfacial  resistance  (Q  cm2) 

Rohm •'  ohmic  resistance  (Q  cm2) 

a:  degree  of  distribution 

Sbet •'  surface  area  measured  by  BET  method 

TEA+:  tetraethylammonium  cation 

TLM:  transmission-line  model 

TLM-PSD:  transmission-line  model  with  pore  size  distribution 
top:  operational  time  (s) 

V0hm'-  ohmic  voltage  drop  (V) 

V totai:  total  pore  volume  (cm3  g-1) 
w:  angular  frequency  (rad  s-1) 

Z;  impedance 

Z interface •'  interfacial  impedance 
Zintra-pore'  impedance  of  intra-particle  pores 


